www.gusucode.com > 循环自相关函数工具箱源码程序 > matlab代做 修改 程序循环自相关函数工具箱/cyclostationary_toolbox/ex3.m

    clear;
clc;
t=0:1/600:1023/600;
A=1;B=1;fb=7;fa=100;
x=A*(1+B*cos(2*pi*fb*t)).*cos(2*pi*fa*t);
figure(1)
plot(t,x);
axis([0 1 -2 2]);
%alpha=2*pi/300;
T=300;
max_tau=0;
R=cyclic_autocorrelation_fast(x,T,max_tau);
y=R.*conj(R);
% ff=600*(0:1023)/512;
ff=0:1:300-1;
figure(2)
% plot(ff(1:150)/(2*pi),y(1:150));
plot(ff,y);